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Abstract : Modern wind turbine technology has been a great 
improvement over the past few decades , leading to large 
scale wind power penetration. The increasing penetration of 
wind power resulted in emphasizing the importance of 
reliable and stable operation of power systems , especially in 
an isolated small-scale power system. This paper describes 
the impact of full-load converter interfaced wind generators 
on power system oscillations. Eigenvalue analysis and 
participation factors are carried out under variant wind 
power penetration levels to investigate the small signal 
stability of the isolated small-scale power system in Penghu 
Island. 
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I. Introduction 

Renewable generation resources have become an increasingly 
popular alternative to conventional fossil fuel based generation 
over the past few decades. Among them wind power is the 
most rapidly developing electricity generation source, and the 
installed wind power capacity has been reached around 396 
GW worldwide by the end of 2015 [1]. As wind energy is 
increasingly integrated into power systems, the reliability and 
stability of wind power generation system has been getting 
more attention from researchers owing to the fact that the 
operational conditions and characteristics of wind generators 
are very different from those of traditional synchronous 
generators [2]. 

Many special intensive researches for the small signal stability 
of wind farms connected to power systems have been 
proposed in the last decade, which devoted to modelling and 
analysis of wind turbines with variant types of generators. In 
[3], models of the three most important actual wind turbine 
concepts: a constant speed wind turbine with squirrel cage 
induction generator (SCIG), a variable speed wind turbine 
with doubly fed induction generator (DFIG), and a variable 
speed wind turbine with a direct drive synchronous generator 
(DDSG), were described and a method was presented to 
calculate the initial conditions of dynamic models of the three 
concepts. The modelling and small signal analysis of a grid 
connected DFIG was presented in [4], and eigenvalue and 
participation factors analysis of the linearized models were 
carried out. Results show that the stator current, turbine speed, 
twist angle, generator speed and rotor flux q-component are 
the main responsible for oscillating modes. In [5], a damping 
controller on the different modes of operation for the DFIG 


based wind generation system was implemented, and an 
optimal algorithm for coordinated tuning of the damping 
controller to enhance the damping of the oscillatory modes 
was presented. The results show that the damping controller 
gives effective damping out the low-frequency oscillations. In 
[6, 7], the small signal analysis considering converter and 
pitch controllers was studied, and a model simplification of 
the DFIG wind turbine was presented to promote dynamic 
performance. An optimal controller design using 

multi-objective differential evolution (DE) was presented in 
[8]. Results show that using the coordinated tuning of the 
controllers’ parameters of the rotor and grid side converters 
can effectively improve the damping performance and the 
stability margin of the DFIG wind turbine. In [9, 10], the 
impact of full-load converter interfaced wind turbines on 
power system small signal stability was studied. The results 
have shown that there is a general decoupling between the grid 
dynamics and the wind turbine mechanical system by the 
full-load converter, that is, it implies that the system 
oscillatory modes of the wind turbine mechanical system have 
less effective than those of the synchronous generators’ 
mechanical system. The eigenvalue analysis based on models 
considering the controllers’ parameters of a wind turbine with 
direct drive PMSG was evaluated in [11, 12]. Mode shape and 
participation factors analysis reveals the origin of different 
modes and controller’s parameters can be chosen properly. 

In this paper, the small signal stability analysis of an isolated 
small-scale power system with full-load converter interfaced 
wind turbines is presented to investigate the impact of 
different wind turbine penetration on power system behavior. 


II. Small Signal Stability 

Power system dynamics can be described by a set of nonlinear 
differential equations together with a set of algebraic equations 
(DAE) of the following form [13]: 
ix = f{x,y,u) 


where* is the vector of the state variables, y is the vector of the 
algebraic variables, and u is the vector of system parameters, 
e.g., the load level of the whole system. In the case of small 
signal stability analysis, (1) is linearized around 
anequilibriumpoint to yield the linearized DAE model as: 


Ax 


0 



f x (u) f y (u ) 
g x (u) g y (u ) 


A* 

Ay 


' ^total 


Ax 
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( 2 ) 


wher eA tota i is the total system Jacobian matrix, Adenotes an 
incremental change in steady state value. By eliminating the 
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algebraic variables y in (2), the equivalent system state matrix 
can be formulated asA = f x (u) - f y (u) g y l (u) g x (u). The small 
signal stability of the nonlinear system (1) in the neighborhood 
of the equilibrium point can be analyzed by inspecting the 
eigenvalues of the linearized system (2); e.g., all of the 
eigenvalues of the system state matrix A should have negative 
real parts to maintain the stability. Accordingly, eigenvalues 
close to the imaginary axis of the complex plane are the 
critical eigenvalues of interest to analyze the stability of power 
system oscillations. 

Equation (3) describes that the critical eigenvalues and 
corresponding eigenvectors construct a dominant invariant 
subspace with respect to the entire eigenspace of system state 
matrix. 

A(u)V(u) = V(u)A(u) 

U(u)A(u) = A(u)U(u) ( ' J> 

where the column vectors, known as right eigenvectors, in V(u) 
and the column vectors, known as left eigenvectors, in U(u) 
are the basis in the invariant subspace, and the eigenvalues of 
A(u) are a fraction of the corresponding eigenvalues in the full 
matrix A (u). 

The participation matrix P combines the right and left 
eigenvectors as a measure of the relation between the states 
and the modes [13]. 


P = \.P\ Pi ■••Pj ( 4 ) 

with 
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v u u n 
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= 
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wherep,is the participation factor of mode i. For the r 
eigenvalue (A,-= ai±fli), the frequency of oscillation in Hz is 
given by (5), and the damping ratio is given by (6), which 
determines the rate of decay of the amplitude of the 
oscillation. 


c t = 


-a: 


a /* 2 ! 2 + A 2 


(5) 
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A 
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III. Description of the Penghu Power System 

Penghu archipelago, covering an area of 127 square kilometers, 
is located in the Taiwan Strait between mainland China and 
Taiwan Island. The Penghu system shown in Fig. 1 is a typical 
weak and isolated power system in which the main power 
generation located at Chien-shan power plant consists of 12 
diesel generators, bring total installed capacity up to 129.77 
MW. Three wind farms, namely Chung-tun WF no.l, 
Chung -tun WF no. 2 and Husi WF, are integrated into the Husi 
substation and have total installed capacity of 10.2 MW. The 
above mentioned power generation provides electricity 
consumption by feeders connected to Magong and Husi 


substations. The demand loads in Penghu system vary 
depending on weather conditions and in general can be 
formulated as four operating scenarios, i.e., summer peak load 
(SP), summer light load (SL), winter peak load (WP) and 
winter light load (WL) for research purpose. For small signal 
analysis, Table 1 shows the models of the diesel generators, 
exciters and speed governors, and the controllers used in the 
three wind farms. Table 2 shows the total state variables and 
demand loads of the four operating scenarios, and 
corresponding wind power generation for the simulation 
conditions of this study. Comparing to the diesel generators, 
the wind generators are expected to have as much generation 
as possible owing to their lower generation costs. Therefore, 
all wind generators operate under the assumed rated power 
output so as to have a high penetration of wind power in 
Penghu power system, especially at light load conditions. 


Magong substation Husi substation 



Chien-shan power plant 


Fig. 1 : Single-line diagram of the Penghu power system 


Models 

Generator 

Exciter 

Speed 

governor 

genenrators #1 -#4 

GENSAL 

IEEEX2 

DEGOV1 

genenrators #5 -#12 

GENSAL 

ESAC8B 

DEGOV1 

wind farms 

full-load converter interfaced wind turbine 

WF controllers 

PQ, current and over frequency power 
reduction controllers 


Table 1: Models of diesel generators and controllers of wind 
farms 


Scenarios 

Diesel 
generato 
r in 

service 

Demand loads 

Wind power (MW) 
/ Number of states 

P 

(MW) 

Q 

(MVAr) 

Wind 
farm in 
service 

Wind 
farm 
out of 
servic 
e 

SP 

#1, #2, 
#5 - #10 

74.807 

7.849 

10.2/340 

0/130 

SL 

#1, #2, 
#5- #7 

31.967 

-4.458 

10.2/292 

0/82 

WP 

#1, #2, 
#5- #9 

62.291 

6.654 

10.2/324 

0/114 

WL 

#1, #2, 
#5- #7 

31.444 

1.371 

10.2/292 

0/82 


Table 2: Number of system state variables, demand loads, 
operational diesel generators, and wind power generation 
under variant scenarios 
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IV. Simulation Results 

The above-mentioned operating scenarios with a varying 
penetration of wind power were investigated and compared to 
the cases with only synchronous generators in operation. 
Eigenvalue analysis under different wind penetration levels 
was implemented to investigate the dynamic behavior in the 
Penghu power system. The diesel generator is modeled with a 
six order state equation, state variables v=[psie, psiD, psiX, 
psiQ, speed, phi]. The exciters ofthe generators #l-#4 and 
generators #5-# 12 are modeled with a six and a five order state 
equations, state variables x=[x a , x h , x e , Xf\, Xf 2 , x r \ and x=[x a , x den 
x e , x h x r ], respectively. The speed governor is modeled with a 
six state equation, state variables x=[xi, x 2 , Jt 4 , v 5 , x 6 ]. The wind 
turbine with full-load converter is modeled with a fifteen order 
state equation, including the states of the PQ, current and over 
frequency power reduction controllers. The detailed 
representation of the state variables can be found in [14]. 
Simulations were realized using DIgSILENTPowerFactory 
15.1 software [15]. 


4.1 Eigenvalues 

Small signal stability analysis identifies the causes of poorly 
damped or unstable power system oscillations and provides 
effective means for oscillation damping control. Eigenvalue 
analysis, based on the computation of the critical eigenvalues 
and corresponding eigenvectors of the system state matrix, 
gives this kind of information for the analysis and control of 
the small signal stability of the system. In this study, the 
eigenvalue analysis has been carried out for the Cases SP, SL, 
WP, and WLinthe Penghu power system, and all eigenvalues 
have negative real parts, which indicated that the Penghu 
power system is inherently dynamically stable under different 
scenarios. The plots of obtained eigenvaluesare shown in Fig. 
2. Comparing Fig. 2-(a) with 2-(b)of each study case shows 
that eigenvalues associated with the dynamics of full-load 
converter wind turbines areactually fast decay modes, i.e. 
eigenvalues with larger negative real parts. It means that the 
dynamics of full-load converter wind turbines has relatively 
limited contribution to the low frequency oscillations. 

The critical eigenvalues associated with the dominant 
oscillatory modes, also known as low frequency oscillatory 
modes, for the Cases SP, SL, WP, and WLis summarized in 
Tables 3-6. Results show that all oscillation modes have 
damping ratios larger than 13%, and only five oscillating 
modes with damping ratios between 13% and 20% in each 
study case. It means that the Penghu power system is a well 
damped system, whether the wind farms are in service or not. 
The observation can be drawn that the Penghu power system 
operating on light load has larger damping ratio than that 
operating on peak load, and the presence of the full-load 
converter wind turbines has less impact on the damping of low 
frequency oscillations. 


Fig. 2: Eigenvalues with/without wind power penetration 
under different scenarios 
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Scenario 

Eigenvalue 

Oscillation 

frequency, 

Hz 

Damping 

ratio 

Summer peak 
loadwith wind 
farm 

penetration 

-2.86359± 

j20.32880 

3.23543 

0.13949 

-3.96518± 

j28.02058 

4.45961 

0.14011 

-3.77347± 
j26. 64261 

4.24030 

0.14023 

-7.10147± 
j46. 98080 

7.47727 

0.14946 

-0.09439± 
j0.3388 1 

0.05392 

0.26836 

Summer peak 
loadwithout 
wind farm 
penetration 

-2.84443± 
j20. 36390 

3.24101 

0.13834 

-3.92909± 
j27. 39341 

4.39163 

0.14097 

-7.10153+ 
j49. 98051 

7.47718 

0.14946 

-4.17473+ 
j24. 90488 

3.96374 

0.16532 

-0.09539+ 
jO. 33791 

0.05378 

0.27168 


Table 3: Dominant oscillatory modes for Case SP 


Scenario 

Eigenvalue 

Oscillation 

frequency, 

Hz 

Damping 

ratio 

Summer light 
loadwith wind 
farm 

penetration 

-3.06665+ 
j20. 36862 

3.24177 

0.14888 

-7.10204± 
j46. 97479 

7.47627 

0.14949 

-3.95243+ 
j25. 26231 

4.02062 

0.15458 

-4.29631+ 
j27. 33522 

4.3505 

0.15527 

-0.09909 ± 
jO. 35009 

0.05572 

0.27232 

Summer light 
loadwithout 
wind farm 
penetration 

-3.04624+ 
j 20.40272 

3.24720 

0.14767 

-7.10224± 
j46. 97429 

7.47619 

0.14950 

-4.26421+ 

j26.6704 

4.24473 

0.15788 

_4.54704± 
j23. 38833 

3.72237 

0.19084 

-0.10195+ 

j0.35501 

0.05650 

0.27601 


Table 4: Dominant oscillatory modes for Case SL 


Scenario 

Eigenvalue 

Oscillation 

frequency, 

Hz 

Damping 

ratio 

Winter peak 
loadwith wind 
farm 

penetration 

-2.88423+ 
j20. 37619 

3.24297 

0.14015 

-3.79408+ 
j26. 26862 

4.18078 

0.14295 

-4.02150+ 
j27. 84183 

4.43117 

0.14296 

-7.10151+ 
j46. 98026 

7.47714 

0.14946 

-0.09394± 
jO. 34369 

0.05470 

0.26366 

Winter peak 
loadwithout 
wind farm 
penetration 

-2.86366± 

j20.41291 

3.24882 

0.13893 

-3.98314+ 
j27. 37512 

4.35689 

0.14399 

-7.10159+ 
j46. 97990 

7.47708 

0.14946 

-4.23051+ 
j24. 52102 

3.90264 

0.17001 

-0.09508+ 
jO. 34301 

0.05459 

0.26713 


Table 5 : Dominant oscillatory modes for Case WP 


Scenario 

Eigenvalue 

Oscillation 

frequency, 

Hz 

Damping 

ratio 

Winter light 
load 

with wind 
farm 

penetration 

-2.97120+ 
j20. 52730 

3.26721 

0.14324 

-4.13310+ 
j 27. 84249 

4.43127 

0.14684 

-3.81696+ 
j25. 65823 

4.08364 

0.14714 

-7.10202± 
j46. 97406 

7.47616 

0.14949 

-0.10771+ 
jO. 35992 

0.05728 

0.28670 

Winter light 
load 

without wind 
farm 

penetration 

-2.96103± 
j20. 55013 

3.27065 

0.14262 

-7.10217± 
j46. 97371 

7.47610 

0.14950 

-4.12252+ 
j27. 14895 

4.32089 

0.15013 

-4.35863± 
j23. 88378 

3.80122 

0.17953 

-0.10886± 
jO. 36414 

0.05796 

0.28642 


Table 6: Dominant oscillatory modes for Case WL 
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4.2 Participation factors 

Low frequency oscillatory modes, usually including inter-area 
and local-area modes, can be identified by the concepts of 
frequency of oscillation, mode shapes and participation factors. 
The participation factors of a mode indicate the contribution of 
the system state variables to the mode. According to the 
equation (4), participation factors of a critical oscillation mode 
were calculated with the critical eigenvalues and 
corresponding left and right eigenvectors of the system state 
matrix, which were based on the four scenarios as shown in 
Table 2. In each scenario, the Penghu power system connected 
with and without wind farms were considered. The results 
have been shown in Figs 3-6 and Table 7. 

Figs 3-6 show that the participation factors associated with the 
most critical eigenvaluehave their magnitudes larger than 0.1, 
and are clustered within 30 degrees in each groupforeach 
study case. Results indicate that the main contribution of the 
oscillations is the state variables associated with the dynamics 
of the synchronous generators in Chien-shan power plant. It 
presents a local electromechanical mode, and is mainly 
associated with generator speed (speed), rotor position angle 
(phi), and flux in Q-damper winding (psiQ) transients. 

In Table 7, comparing the penetration factors with and without 
wind farms penetration for each study case shows that all of 
the state variables associated with the full -load converter 
interfaced wind turbines do not contribute to the local 
electromechanical mode, it means that the configuration in the 
full-load converter wind turbine has the ability to decouple the 
wind generator dynamics form the gird dynamics. 



Fig. 3: Comparison of the participation factors of the critical 
eigenvalue for case SP 



Fig. 4: Comparison of the participation factors of the critical 
eigenvalue for case SL 
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Fig. 5: Comparison of the participation factors of the critical 
eigenvalue for case WP 



Fig. 6: Comparison of the participation factors of the critical 
eigenvalue for case WL 


V. Conclusion 

The small signal stability of the wind farm integrated power 
system was examined. Based on the eigenvalue and 
participation factor computation, the simulation with various 
wind penetration levels and wind turbine with full-load 
converter was conducted to examine the low 
frequencyoscillations of the Penghu power system. The 
simulation results show thatthe eigenvalues associated with 
the dynamics of the synchronous generators are the mainly 
local electromechanical modes, whereas those with full-load 
converter wind turbines are actually fast decay modes and thus 
have less impact on the damping of low frequency 
oscillations. 
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Partic . . 

c+ . . .. Participat 

A State ipatio . 

Cluste T t , ion 

Vanabl n . . 

r r Angle 

e Magn (del) 

itude , g) 

^ . . Partic 

0 Particip . 

State . ipatio 

Clust T 7 ... ation 

Vanabl __ . n 

er Magmt A 

e . Angle 

Ude /A \ 

(deg) 

Case SP with WFs: X = -2.86359 
±j20. 32880 

Case SP without WFs: X = 
-2.84443 ±j20.36390 

0 331 

G5_ phi 134.471 

1 1 

G6_phi 0485 1 32.3 85 

1 G10; 0.439 99.486 

speed 1 

G7; phi 0.996 167.181 
9 

G9; 0.596 168.115 

speed 0 

G6; phi 0.238 151.461 
4 

0 996 

2 G9; phi ^ 167.051 

G10; 0419 -H3.044 

phi 8 

G7; phi °' 276 -100.118 

3 ooss _ 99141 

psiQ 0 

G8 ’ 0414 -131.273 

speed 7 

G10; 0.437 -93.331 

phi 7 

G8; 0.249 -121.47 

speed 6 4 

G7; 0.158 -99.055 

psiQ 9 

G1 °’ °' 421 79.905 

speed 1 

G7; °' 277 92.831 

speed 3 

4 G8; phi 0.248 45.707 

8 

G7; 1 0 

speed 

G9; phi 0.594 -24.702 
5 2 

G6; 0.239 -15.719 

speed 1 

5 G8; phi 0414 35.778 

6 ° 9 ’ 1 0 
speed 

G5; 0.332 _ 32 5gl 

speed 1 

G6 ’ 0485 -34.666 

speed 9 

Case SL with WFs: X = -3.06665 
±j20. 36862 

Case SL without WFs: X = 
-3.04624 ±j20.40272 

G5; phi 0.365 101.679 

1 8 

G7; 0.225 87.3925 

speed 2 

1 G6; 0.107 122.923 

psiQ 4 

2 G7; phi 0.995 166.149 

0 

2 G6; phi 0.995 166.011 

0 

G5 ; phi 0.975 -111.48 

9 9 

G6; 0.522 -137.89 

speed 0 3 

G7; 0.205 -99.182 

psiQ 8 

G7; phi 0.224 -106.595 

3 1 

G6; 0.206 -99.269 

psiQ 5 

4 G6; 1 0 

speed 

4 G5; 0.980 82.361 

speed 7 

5 G5; 0.367 -64.332 

speed 7 

G7; 1 0 

speed 

5 G6; phi 0.519 28.256 

52118 
7 


Case WP with WFs: 2 = 
-2.88423 ±j20.37619 

Case WP without WFs: X = 
-2.86366 ±j20.41291 

G5; phi 0.996 166.979 
5 

G9; phi 0.951 119.184 
3 

] G7; 0.659 -171.884 

speed 6 

G6; phi 0.235 138.202 
3 

G9; 0.160 -146.878 

psiQ 1 

1 G8; 0.158 109.321 

speed 6 

2 G9; phi 0.996 167.119 

6 

G6; phi 0.712 -157.79 

3 ° ° 

G7; 0.206 -153.96 

speed 7 1 

G9; 0.167 -98.992 

psiQ 5 

, G8; phi 0.158 -83.559 

4 i 

G6; 0.119 -63.901 

psiQ 6 

2 G5; 0.167 -99.082 

psiQ 7 

3 G7; 0.110 89.032 

psiQ 6 

G5; 1 0 

speed 

G9; 0.954 -47.795 

speed 5 

G7; phi 0.657 -4.904 

3 

G6; 0.236 -28.777 

speed 1 

G9; 1 0 

speed 

G6; 0.714 35.090 

speed 4 

G7 ; phi 0.206 13.157 

0 

Case WL with WFs: X = 
-2.97120 ±j20.52730 

Case WL without WFs: X = 
-2.96103 ±j20.55013 

G5; phi 0.994 166.648 
9 

G6; phi 0.248 166.562 
1 8 

G7; phi 0.248 166.562 
8 

G6; phi 0.755 124.189 

1 8 

G5; 0.185 95.037 

speed 0 

2 G7; phi 0.995 166.728 

3 

2 G5; 0.205 -98.820 

psiQ 1 

3 G6; 0.154 -141.34 

psiQ 3 9 

G5; 1 0 

speed 

G6; 0.249 -0.0322 

speed 8 

G7; 0.249 -0.0322 

speed 8 

G7; 0.203 -98.810 

psiQ 2 

G5 ; phi 0.184 -98.233 
2 

5 G7; 10 

speed 


Table 7: Participation factors and associated state variables of 
the critical eigenvalue for each Case: SP, SL, WP, and WL 
under with/without wind farms penetration 
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